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closely resembled observed oceanic variations. Atmospheric boundar: 
layer model reformulation is required to alleviate the cloud 
thickness/shortwave radiation problem. 
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ABSTEACT 


A Coupled, océanic-atmosfpheric boundary layer model 
which provides single-station prediction capabiiity is eval- 
uated relative to boundary layer observations. The model is 
initialized ana verified using data obtained during the 1983 
Mixed Layer Dynamics Experiment (MILDEX). Model prediction 
of inversion height, lifting condensation level, air and sea 
temperatures, specific humidity and mixed layer depth are 
compared with observations. A Significant model shortcoming 
is the over-prediction of cloud thickness. Corseguently, 
Shortwave radiation at the ocean surface is toc low and tne 
predicted ocean mixed layer depths are rot freaiistic. 
Oceanic predictions resulting from independent svecification 
of surface radiation more closely resembled observed oceanic 
variations. Atmospheric boundary layer model reformulation 
is required to alleviate the cloud thickness/short wave ta aie 


ation problen. 
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fT. INTRODUCTION 

All aspects of naval operations are affected by the 
adjacent ocean/atmosphere envircnment. They perrortanecesen 
weapons systems operating on the principles of electro- 
optical (EO) and electromagnetic (EM) physics are dependent 
upon the existing atmospheric conditions in the marine 
atmospheric boundary layer (MABL). Those systems whose 
operation iS yoverned by water borne acoustic paths are 
dependent upon the environmental conditions in the ocean 
boundary layer. Aside from the exceptionaliv severe weather 
phenomena so commonly associated With environmental 
dynamics, the complex weapons systems in use througnout the 
fleet are affected by routine environmental conditions and 
subtle changes therein as well as by the renovned extrenes. 

Atmospheric dynamics wili alter the vertical profiles of 
temperature, humidity and pressure which are the variabies 
affecting the refractivity of the atmosphere. The vertical 
gradient of refractivity is the critical factor in the pres- 
ence of phenomena such as ducting, Subrefraction, super- 
refraction, frequency dependencies and their impact upon tke 
performance of electro-optical and electromaynetic svstems. 

Oceanic dynamics manifest as vertical variations in 
temperature, pressure and salinity affect the sound speed 
Phoz tle. As the most critical factor in determining the 
propagation of acoustic energy, variations in the sound 
speed profile as a function of depth will result in 
phenomena Such as surface ducting and convergence zones and 
Will determine cutoff Frequencies and optigaal source/ 
receiver geometries. 

Atmospheric and oceanic bcundary layer processes are 
nighly interdependent. The ability to prediete their 1 ntem. 


action is clearly essential to the tacticai commander. 
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The model evaluated consists of a one dimensionai atnos- 
pheric model (Davidscn et al., 1984) and a one dimensional 
oceanic Eoundary layer model (Garwood, 1977). Tne models 
were coupled (O'Laugniin, 1982) and demonstrated to realis- 
tically simulate boundary layer idlynamics (Hervey, eons). 
The goal of this study is to assess the ability of a coupled 
ocean-atmosphere boundary layer model to accurately take 
such predictions by conducting a boundary ilaver simulation 
Usama cecamec and atmospheric data obtained Guring the 1983 


Mixed Layer Dynamics Experiment (MILDEX). 
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II. BACKGROUND 


A. MOTIVATION 


The evolution of the NPS coupled atmoSpheric and oceanic 
boundary layer model has reached the developmental stage of 
being used to predict atmospheric and oceanic dynamics using 
actual data sets for initialization and comparison with 
Od ao meo wr. Prior efforts along this model development 
path include an examination of the coupled model's anpility 
to forecast uSing real data for initiaiization and coapar- 
ison or results in order to determitie the adequacy of the 
coupling (C'Latghiriye i Js2) This study concluded that the 
coupled model did indeed yield fairly accurate scesuits; 
however, the data sets employed (CEWCOM 1976 and CEWCOH 


1978) were inadequate. The inadequacy was a resuit of the 
frequent transit oof the data gathering vessel, RV ACANTA 
over Significant distances through ainnomogeneous water 
masses. 


A follow-on study examined the adeguacy of the coupled 
boundary layer model reiative to the performance of the 
uncoupled atmospheric and oceanic models. Tue results indi- 
cated that the benefits of coupling to the MABL were mnino7z; 
however, the predictive capability of the OBL was greatiy 
enhanced (Hervey, 1983). 

Both studies recommended evaluation with data sets where 
effects of aivection and vessel movement are minimized. 
such a data set was obtained during the 1983 Mixed Layer 
Dyhamics Experiment (MWILDEX) (Geernaert et al., 1984), 
conducted from October 24 through November ets oo. Data 
were obtained 3 to 5 November, 1982 Tne ie eiyeet a1 ty enue 
and 126 W ky RV ACGANTA (began Meteorological data 
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Figure 2.1 Positions Of RV ACANEA: 03-05 Nov 1983 


consisted of atmospheric pressure, windspeed, temperature, 
humidity and radiation in the atmospheric surface layer. 
Profiles of atmospheric wind speed and direction, tempera- 
ture and humidity using a navaid radiosonde system were also 
obtained. Additionally, observations of sea state and 


cloudiness were recorded hourly. 
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Cceanographic data were obtained with a Neil 3Brown 
Instrument System (NBIS) CTD. The device recorded values of 
temperature, conductivity, depehgaid dissolvecdyom geq. Two 
other research vessels, WYACOHMA and FLIP acquired data in the 
Vicinity sof eeu Naan FLIP was within 2 kilometers o£ the 


ACANIA from October 26 through November 9. 


B. BOUNDARY LAYER CHARACTERISTICS 


The Marine Atmospheric Boundary Layer ({MNABL) extends 
from the air-sea interface through a capping inversion tvpi- 
calidy 190 to 100 0n signe The Oceanic boundary Laver (O3L) 
extends from the interface to the thermocline, 2 to 50m 
deep. The MABL is comprised of a cool, moist, turbulent 
well mixed layer (relative to the upper atmosphere) in which 


equivalent potential temperature and specific humidity are 


basically constant with height. At the inversion levei a 
Significant change (jump) in the vertical profiles of these 
variables occurs. The air above the inversion iS warmec 


than that of the mixed layer and is characterized Ey a lapse 
rate of increasing temperature as a function of height. 
Similariy, the air of the free atmosphere is drier than that 
of the MABL with the moisture content decreasing as a func- 
tion of heignt. 

The value of the well mixed quantities change over time 
due to the turbulent fluxes of heat and moisture at the 
air-sea interface and at the inversion. Turbulent Kinetic 
energy provides the necessary forcing of thase processes. 
The turbulent kinetic energy budget provides the foundation 
upon which most of the modelling is based. The weli mixed 
nature of the variabies in the MABL and OBL imply that heat 
and moisture are conserved during the mixing process and 
that the variation with height oF the vertical turpulent 


rluxes of momentum, humidity and neat is linear. Assuming 
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advection is negligible, these properties permit the 
modeiing of the variation of the variables based solely uvon 
surface and inversion layer height fluxes. 

Turbdulence at the inversion causes the entrainment of 
WacLm, dry air into the MABL and growtn of the Layer, large 
fsyuepeic) scale atmospheric forcing contributes, to subpsi- 
dence. The time rate of change of the height of the MABL is 
aesesultevof these effects. 

An important atmospheric feature is the presence of 
elouds. If the lifting condensation ievel (LCL) is lower 
than the height of the inversion, clouds will be present. A 
Sum@dowcoetayer in the MABL will vary the amount or solar, 


Shortwave and return long wave radiation in the boundary 


layer. Such variations will greatly affect the dynamics of 
the NABL and OBL (Hervey, 1983). 
The Cceanic Boundary Layer (OBL) 1s a turbulent region 


of the upper ocean bounded by a dynamically stable water 
mass at its base and by the air-sea interface above. Tae 
layer is characterized by constant velocity and density 
profiles arising from the high degree of turbulence within 
the layer anda relatively slight amount of turbulence in 
the underlying stable watermass. Jomentum and buoyancy 
fluxes through the surface provide the energy which resuits 
in the boundary layer turbulence. 

Density is a function of temperature and salinity, both 
of which are assumed to be homogeneous in the boundary 
faye « Since the model will forecast for relatively short 
periods of time (less than 72 hours), the effects of 
Salinity are assumed to Fe negligible (Miller, 1376). As in 
the case of the MABL, the vertical homogeneity of the OBL 
allows for the modeliing of dynamic processes through the 
use of only the surface and entrainment fluxes. The turodu- 
ient kinetic energy budget again serves as the foundation 


for nodeling. 


ie 


C. #WIXED ATHOSPHERIC BOUNDAKY LAYER MODEL (MABL) 


The atmospheric variables required as input are short- 
wave radiation, which depends on Julian date, time cf day 
and latitude, inversion height, mixed layer and jump vaiues 
for eguivalent potential temperature and specific humidity. 
Also the above layer lapse rates for specific humidity and 
potential temperature, subsidence rate, wind speed and 
direction within the mixed layer and sea surface temperature 
are required. In the uncoupied model sea surface tempera- 
ture is constant. However, the coupled model uses the vari- 
ations in sea surface temperature as predicted by the OBL 
model aS an input. 

Three methods are availakle for the computaticn of 
sutsidence (large scale vertical velocity}. The methods are 
the kinematic, the adiabatic and the integration of the 
moisture budget equation (Q method). The 9 method is the 
most accurate approach given Single station data, according 
to Gleason (1982). It was used for subsidence calculation 
for Ghis. study 

The MABL predicts the evolution of the LCL, inversion 
height, mixed layer values of equivalent potential tempera- 
ture and specific humidity, their respective jump values, 
entrainment rate, short and Vong wave Gadiacle a eee le en 
velocity (U*), and the scaling parameters for temperature 
and humidity (T* and Q*) for a fpreScribed period. The model 
uses a 30 minute time step. The computational steps are 
depleted inp dea 2. ce 

The time evolution of the conservative variables are 
predicted using integrated (with respect to height) sate 
eguations (Tennekes and Dreidonks, 1981). These equations 


are: 
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Pigure 2.2 MABL Schematic 


Ve 


h (DX/Dt)= (w' x") (W'x") +source (26. In) 
h (DAx/Dt)=h Vx (oh/St)~ (w'x") + (w'x') ~source (252) 


where the subscripts zero and h refer to surface and inver- 
sion height values, respectively. Gamma (tT) 1s the lapse 
rate above the inversion height and the sources term is egual 
to - (Fan ~Fno)/?Cp for X=©@ and egual to zero for X=q. The 
Subscripts h and zero refer to inversion height and 
surface, respectively, Fn is the net radiative flux. The 
quantity (sh/dt) is the difference between tne time rate or 


change of the inversion height (dh/dt) and subsidence rate. 


Tne system of eguations is closed through the use of 


Stage and Businger's (1981) entrainment rate parameteriza- 
Galen: This is also used in the determination of the 
temporal evolution of the inversion heignt. To obtain the 


entrainment rate, Stage and BusSinger made a clostre assump- 
tion that the dissipation rate of turbulent kinetic energy 


is a fraction ({1-A) of the production Hace, 7. 


D= {1-A)P (2.3) 


The quantity 'A* in essence is the critical flux Richardsor 
number (Ricr) which separates turbulent from nonturbulent 
boundary layers. It has been determined through the 
modeling of dry convective layers that A=0.2. This value is 
in close agreement to that calculated by Businger (1973) and 
lies in the middle of the range of critical flux Richardson 
numbers, 0.15 to 0.25, which were determined by Arya (1972). 

Earlier modeling efforts indicated that only 1to 4 
percent of the TKE produced in the boundary layer was used 
ror entrainment. The 0.2 factor however indicates that 20 


percent of the TKE is thus utilized. 
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Bulk aerodynamic formulas are usSeji in the determination 


of the surface fluxes of momentum, heat, and ndisture: 


*=ca *D (2.4) 
T*=Ce 4e-@) (2.5) 
g*=Ce"(g-qo) (Ach) 


The fluxes are represented by: 


re) -u+? ear 
(-T'w!)=U*T* (223) 
(-g'w")=U#q+ (2.9) 


Cd and Ce are the stability dependent drag and exchange 
coefficients. The zero subscript refers to surface values, 
overbars represent average quantities and primes (') repre- 
sent fluctuations from mean values. 

Long and shortwave radiation fluxes ace calculated inde- 
pendentiy. Modeling of these quantities is extremely diffi- 
Cie. Effects of absorbing aerosols and gases, and of the 
distribution of cloud droplet size, are nearly impossible to 
predict . Long wave radiation modeling is based upon empir- 
ical relations between the average cloud liguid water 
content and cloud emissivity. Cloud top net long wave radi- 
ative flux L is determined from the Stefan-S5o0ltznan iaw and 


cloud top temperature: 
L= oO (Th *- Th) (2. 10) 
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The parameter § is Stefan's constant and the parameter €cis 


determined by: 


Ec=1-exp (-a¥) (2211) 


where & is the cloud emissivity, Wis the cloud liguid water 
content, and tat is a constant equal to .158 n2g! (Slinger 
et al., 1987 Nee Cloud base long wave radiative flux is 


represented by: 


Ee Ce Soe eey (25512) 


where Ts is the sea surface temperature and Tc is the cloud 
base temperature. Radiative flux divergence inthe air 


between the sea surface and cloud base is neglected so that 


Li =Lage Therefore, the net long wave flux at the surface is 
given by: 
F (long) =0(Ts4-@T4= (1-&) Te) (2273) 


where T is the average cloud temperature. 

For the cloud free condition, the net lonj wave flux is 
calculated at z=h and z=0 by integrating tne flux emissivity 
profile as developed by Fleagle and Businger (1930). Quite 


Simply, the net long wave radiative flux is given by: 


F (long) =Fu-Fd (2.14) 


where Fu and Fd represent the net upward and downward radi- 
ative fluxes, respectively. 

Short wave radiative flux is calculated uSing the 
delta-Eddinston method described by Joseph et al., (1976). 


2,0 


The solar tlux is evaluated over 15 equally spaced {.1um) 
bamaosrangingertrom 0.2 to 1.7 um wavelength. Solar zenith 
angle as determined from previously described input vari- 
ables 1s an essential part of these calculations. Several 
considerations are essential in modeling shortwave radiative 
iiax. Shortwave extinction is primarily a result of scat- 
tering (instead of absorption). Scattering which creates a 
diffuse source of radiation is due to atmospheric particles, 
water dropiets (in the case of clouds) and sea aerosols. 
The MAPL incorporates the direct and diffuse sources of 
shortwave radiation in determining the total surface flux. 
Sea surface rerlection of incident shortwave radiation is 
assigned a value of 10 percent. 

Due to the one dimensionality of the model and several 
assumptions incorporated in the Q method of subsidence 
Getermination, atmospheric advection was assigned a value of 


Z2eLO-. 


D. OCEANIC BOUNDARY LAYER MODEL (OBL) 


The OBL model was developed by Garwood (1977). Pt wee a 
one dimensional, second order bulk model. The nodel employs 
the Navier-Stokes equation of motion with the geostrophic 
component eliminated, the continuity equation (assuming 
incompressible water), the heat equation derived from the 
first law of thermodynamics, a noniinear eyguation of state 
and the conservation of salt equation. 

The dynamics of the entrainment process determine the 
rate of deepening or shallowing (retreat) of the mixed 
layer. The fundamental assumption is that the turbulence of 
the mixed layer provides the energy required to destabilize 
and erode the underlying stable water mass (Garwood, i977). 


AS a result, the turbulent kinetic energy budget provides 
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the basis for the determinaticn of tne entrainment rate. 
The system of equations is closed via the mean turbulent 
field modeling of the vertically integrated equations for 
the individual turbulent kinetic energy components along 
with the bulk buoyancy and momentum equations. 

The OBL model differs from earlier efforts in the 
following ways. First, the amcunt of wind generated turbu- 
lent kinetic energy to be used in mixing is a function of 
the ratio of the maxed layer depth to the Obukhov mixing 
length. Also, viscous dissipation is dépendent on a local 
Rossby number and seperate vertical and horizontal equations 


for turbulent kinetic energy are used (Garwood, 1977). 


Buoyancy conservation employed in the OBL is a generali- 
Zation of the concepts of heat and salt conservation. the 
buoyancy equation is a combination of a linearized equation 
or state, conservation of heat and conservation of salinity 


eduations. Buoyancy is Simply defined by: 


b=(t-0 ) g/ee (2. 15) 


where ®@ is given by: 


© = ko [1-«(8-90) +Q (S-So) } CZ on) 
and where 6,5, and are instantaneous values of temperature, 
Salinity and density, “ and 8 are expansion ve9craaer ents fou 
heat and salt respectively, The zero subscript denotes a 
constant, representative value. 

Parameters required for OBL model initialization include 
Shortwave radiation absorbed from the ocean surface to a 
depth of one meter, the radiation extinction coefficient for 


determining the dewnward heat fiux and thesncriemed 


a 


Richardsen number for a stability adjustment at the base of 
the mixed layer. Surface boundary conditions recuire speci- 
fication of wind speed and direction jlata, cloud cover, sea 
surface temperature, dry bulb and dew point temperatures, 
incident solar radiation and precipitation (PP! and evapora- 
Elon (f) rates. 

Bulk aerodynamic formulae are used to calculate the 
turbulent fluxes of latent heat (Qe) and sensible heat (Ch). 
Shortwave radiation is calculated as previously described 
(delta Eddington method). An empirical eguation (Husby and 
Seckel, 1978) 1s used in the determination of the net back 
Bereatton, Ob. The upward heat flux 1S simply the sun of 


these rluxes: 


Qu=Ce+Ch+Ob C2. lea) 


The net radiation is not simply the difference between the 
downward and upward fluxes since only a fraction of the 
incident radiation penetrates the mixed layer. The depth of 
penetration is a function of turbidity, that is, in coastal 
regions more energy will be absorbed than in the less 
turbid, open ocean regions. The fraction of Shortwave radi- 
ation which is absorbed 1S conceptually regardeaq as a 
GOnmLEIDaton to the upward radiation. Therefore, tke net 


heating at the surface is represented by: 


Onet=Qu+t+REF (Qs) -Os (Ze 1) 
where kf is the fraction of incident radiation absorbed in 
the first meter. 


Given Qnet the surface fluxes of buoyancy and monentun 
can be calculated. MeBeriseDGrntsene turbulent fluxes of 


puoyancy, salinity and temperature can be determined by: 


2S 


b'wi=g (A(T'w")—t(s "wag (22S) 
s'w'=(P-E) So (2520) 
Diy t= One ep (2229 


where So represents a surface Salinity. 

Temperature and Salinity profiles as well as the wind 
driven horizontal current profiles are reguired for OBL 
model initialization. Temperature 1S considered to be the 
dominant factor in density variation when compared to the 
effects of salinity over relatively short time periods such 
aS those encountered in this study (Miller, tees The 
initializing profiles are used in the determination cf mixed 
layer depth. According to tne OBL model, mixed layer depth 
is defined to be the shallowest depth at which the observed 
density value 6t, is .026 t units greater than the density 


at tne surface. 


E. COUPLED BOUNDARY LAYER MODEL 


The processes outlined in the preceding sections have 
been combined to forma coupled boundary layer model. The 
coupling iS achieved by matching the atmospheric and 
oceanic fluxes of momentum, sensible and latent heat and 
radiation at the air-sea interface, O'Laughlin (1982). 

The only oceanic variable which forces the atmospheric 
model is sea surface temperature. Peedback occurs between 
several atmospheric variables and sea surface temperature. 
These include wind stress andthe fluxes outlined previ- 
ously. Conversely, neatly all of the vatmteootemic vdevarice 


influence the evolution of the cceanic variables. 


24 


leeeestedce! tTnitialazation 


Coupled model initialization requires the Julian 
date, local time of day and latitude. Required atmospheric 
vactiables are surface pressure, inversion heignt, LCL, and 
mixed layer values of specific humidity and potential 
temperature. Also reguired are the inversion "jump" values 
and above layer lapse rates of potential temperature and 
Seecific humidity, vertical velocity (subsidence), and 
values for temperature and moisture advection. Tne 10 meter 
wind direction is also specified. The atmospheric nodei 
Simply requires wind speed; however, the wind stress calcu- 
lations needed for the OBL require the direction from waich 
the wind is blowing. These computations lead to the deter- 


Mination of ocean turbulent velocity flux U*? where: 


i - Qaute= 0 wwe (Os 


Oceanic variables required for initialization are 
mixed layer depth, sea surface temperature, temperature jump 
value at the thermocline and the below layer thermal 
gradient. 

Since atmospheric processes occur ona shorter time 
scale than oceanic processes, the MABL model is designed 
with a 30 minute time step whereas the O5FL model uses a 60 
Minute time step. It is assumea that sea surface tempera- 
ture will not change significantly over a 60 minute period. 
As a result, the ocean model is incorporated as a subroutine 
of the atmospheric model and is called every cther time 
step. A schematic of the coupled model is depicted by Fig. 
Ze 
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hee ee PROACH 


The data used in this study are from the Mixed Layer 
Dynamics Experiment (MILDEX), conducted in the Eastern 
Pacific from 24 October to 10 November 1983. Several 
criteria were used in selecting the periods for which a 
sensitivity analysis could be conducted. Tnese included an 


absence of atmospheric fronts and the availability of 


frecuent and accurate radiosonde and CTD data. Frontal 
effects; such as, moisture and heat advection are not 
included in the coupled or uncoupled models. The presence 


£ such activity would therefore impair comparison of the 
one-dimensional model computaticns with observed changes in 
the ocean and the atmosphere in particular. 

The necessity of freguent radiosonde and CTD data is 
critical to the determination of the temporai variation of 
atmospheric and oceanic variables. The radiosonde data had 
to be calibrated using measured temperatures at 20 m, which 
was assumed to be continually within the atmospheric mixed 
layer. The measured temperatures were converted to poten- 
tial temperature, assuming a hydrostatic balance. The 
atmospheric model was then initializei with this potential 
temperature profile. Several data transmission problems 
limited the amount of acceptable radiosonde data. 
Radiosonde launch times coincided with weather satellite 
overhead time. This time resolution was insufficient to 
resolve well any synoptic variability in the atmosphere. 
CCeotwenOarasnegulred no post Calibration. However, it too 
was limited in continuity due te equipment nalfunction. AS 


a result, two periods were selected for this analysis. 


Za 


The first period began at 1600 (Pacific Day ive ees 
November and concluded at 1554 04 November. The second 
period began at 1554 04 November and concludei at 160) 05 


November. The data selection process 15 outlined in Fig. 
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Figure 3.1 Period Selection Criteria 


Initial model values for the atmospheric mixed layer, 
inversion jump and above layer lapse rates for potential 
temperature and specific humidity were obtained from radio- 
sonde data. Sequential atmospheric soundings were used to 
determine subsidence rates, Methods For determining supdsi- 
dence uSing mean values as well as least squares values are 


available. Significant differences were not found for the 
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different approaches; thus, the mean value method was 
employed in model initialization. Sequential radiosonde 
profiles preceding and including the profiie at the begining 
or each experimental feriod were used in subpsidence calcula- 
tikens . The remainder of the atmospheric variables required 
in initialization were tabulatei by Geernaert et al., 
(1984), with the exception of LCL. This value was deter- 
Mined through the use of a Skew T, log P diagram. 

Oceanic variables reguired for initialization were 
determined from temperature profiles from CTD casts. Mixed 
layer depth, temperature jump at the thermocline and below 
layer gradients were thus obtained. The overlapping of 
oceanic data from ACANIA CTD measurements and FLIP provided 


a temporally complete data set. 


B. CASE I: 1600 PDT 03 NOVEMBER-1554 04 NOVEMBER 1983 


The values used to initialize Case T are listel in 
Table i. The uncoupled and coupled model predictions are 
presented in Figs. 3.2 and 3.3 respectively. In terms of 
the evolution of the inversion height and LCL, the unccupled 
model began with a clcud layer (15 m thickness) with a base 
at 266 m. Within two hours, the LCL decreased rapidly (from 
266 to 126 m) as the inversion height rose. The result was 
an extremely thick cloud which attained a maximum thickness 
of 635 mn with a base at 183 m, at 09500. 

The decrease of the LCL was a result of a combina- 
Emon. Of p factors. The specific humidity increased slightly 
after initialization followed by a very gradual decrease 
(-55 g/kg) through 0200. Air temperature decreased eee 
from 1600 through 0300. As a result, the air-sea tempera- 
ture difference increased since the sea surface temperature 


was held constant. The increase in air-Sea tenperature 


ae 


HEIGHT (M) 
O(G/KG) TEMPERATURE 


SPEEOUIZS: 


DEPTHItM) 


1000 


INV(LINE) 


LCL {OASH) 


20 
AIR ee Cee) 


SEA TEMP (DASH) 


SPEC HUM 


WINO SPEED 


-20 Mb xeOs LAY ERMDER iH 


-40 
1600 2000 O0000 0400 A800 1200 1600 


Time 


Figiee 3.2 OUncoupled Model-Case I 


30 


HE 1 GHT (4) 
SPEED(IM/S] O(G/KG) TEMPERATURE 


wee ot ee 
w Oo ~— Nm 


DEP TH(M) 


1000 
INV(LINE) 
S00 
LCL (OASH) 
20 


pieeweniP (LING) 


Sohail (Ono 


SPEC HUM 


WIND SPEED 


ae MIXED LAYER DEPTH 


-50 = 
}600 2000 0000 0400 0800 1200 1600 
Time 


Figure 3,3 Coupled Model-Case I 


S| 


$a 


TABLE I 


Initialization Values-Case I 


Julian Date Sfc Pressure (mb) Sea Sfc Temp (°C) 
307 102224 Vo. 76 

Inversion Ht (m) Mixed Layer Potential Temp (°C} 
281 eo 

Air Temp Jun a6 Temp Lapse Rate (°C/m 
5 CBP Pp (°c) Pp vapes ( ) 


OOS oceans (g/kg) Spec oe Jump (g/kg m) 


Svec Humidity Lps Rate /kgq m Lek. An 

P 76035 y Lp (g/kg m) as ) 

Subsidence (n/s Mixea Laver Depth ({(n 
Sots eee yer Crea 


Ocean Temp Jun oC Ocean Temp Gradient (°C/n 
; Pp p(s) Eaten { ) 


eee A i a I em il a fy A a a a ee a ce 


Nes aes eae eee ees oocege cme no aes = ee a ees seman toeamies ES ee oti Ubi OS a es CES gent, CS a ge) (pan eer tte es Qa (ERO y eee ey | 


Wind: Time Speed (m/s) Dale aed 

1600 om O40 

1802 4.2 319 

2105 8.6 O26 

0017 as Ego 

0417 6.0 276 

0722 Ses 293 

Ug 72 6 a2 028 

Vi 23 6.9 342 

1340 6.1 001 

1541 5.4 340 
difference is reflected 2h a continuous increase in 
MOlSt ure  fiax, Following sunrise, the air temperature 


increased associated with a decrease in specific humidity, a 
less rapid increase of moisture flux and an elevation of the 
Lol 

The inversion height increased throughout the 
prediction period. The time rate of change of this paramn- 
eter is the difference between the rates of entrainment and 
subsidence. Entrainment rate reached a maximum value of 2.3 


cm/s at 2130 and decreased slightly to a value of 1.5 cn/s 


a2 


at the end of the period. Subsidence was held constant at a 
value of -.0032 m/s throughout the period. This value is 
representative for subsidence in this region as noted by 
Davidson et al., (1984). As a result, the evolution of the 
inversion height was dominated by the rate of entrainment. 
The entrainment of warm, dry air into the mixed layer also 
contributed to the gradual decrease of specific humidity and 
the increase in air temperature observed prior to sunrise. 
Sea surface temperature ard mixed layer depth remain 
constant with the uncoupled model execution i.e. no ocean 


variakiiity is permitted. 


2. Coupled Model 


The uncoupled and coupled models predicted nearly 
identical atmospheric results. Clouds were present from the 
beginning of the period. Cloud thickness and the evolution 
of the LCL and inversion height remained the sane as in the 
uncoupled model. Air temperature and specific humidity 
profiles were also unchanged. The coupled model determines 
sea surface temperature hourly. Thererore, the air-sea 
temperature difference, moisture and temperature fiuxes 
would be expected to vary. In this case, the sea surface 
temperature decreased from an initiai value of 18.78 “C to 
18.69°C, a change of only NOG Gaeta saeaesilt /thertime race 
of change of moisture and temperature fluxes for the coupled 
and uncoupled models is nearly identical. Sea surface 
cooling results in convective mixing (negative buoyancy 
Suc cooring throughout the layer, and an increase in 
mixed layer depth for both the atmosphere and ocean. The 
impact of the ocean upon the atmospheric mixed layer 
processes is almost invariant in this case. The major 
influence upon the atmospheric mixed layer is the inversion 


height growth and associated entrainment. 
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C. OBSERVATIONS VS BeGpereoirpua 


Observed specific humidity compares favorably with model 
OUEDUt, FLG eee The time series are nearly identical 
eight hours following initialization. From %2 to 18 hours 
after initialization the observed specific humidity shows an 
increase of approximately 1.0 g/kg. The model predicted a 
continuous decrease throughout the period. As the end of 
the period iS approached, the model and observed specific 


humidity differ by less than .5 g/kg. 


LEGEND 
OBSERVED 
UNCOUPLED 
CUUPLEL 


0 

_ 5 

i 0 

Pe 

o) S 

= 0 

oO 

e 

als 

Tee) 

2 0 

aie) 

£ S 

© 

eae) 0 

ae S 

i 

o 0 

io é 
0 





Ba guee 3c Specific Humidity 
peeked vs Model (Case I) 


The model prediction of air temperature differs by as 


nuch as 2°C from observed values from 2000 to 0400, Fig. 
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Figure 3.5 Air Temperature 
Observed vs Model (Case I) 


Sao < From 1000 to the end of the period, the observations 
and predictions demonstrated an increase and eventually 
attained a constant air temperature, differing by less than 
oe aoe The large difference in air temperature (ca. 2 C) 
which occurs was a result of thick cloud cover and resultant 
shortwave radiation deficit. 

There were vast differences in observed and modeled 
Swemewave bpad@wation, Fig. 3.6. The radiation for a cloud 
free case verifies that there was a cloud cover present. 
This result is deduced from the difference in shortwave 
radiation between the cloud free case and the observed 


values. Glouds were present from 15 to 20 hours after 
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Pigure 3.6 Shortwave Radiation 
Observed vs Model (Case I) 


initialization. Uncoupled and coupled model calculations of 
shortwave radiation are closely correlated. However, they 
are substantially smaller than the observed values. This 
difference 1s a result of the assumption of complete cloud 
cover when the LCL is beneath the inversion height. Horizon 
to horizon cloud cover will result ina fractional amount of 
Shortwave radiation in comparison with a variable amount of 
cloud cover. 

The effect of the over prediction of cloud cover has 
been addressed previously. The cooler” “aire: temperature 
determined by the model and the cooling of the sea surface 
and deepening of the mixed layer by the coupled model are 


two examples. These differences impact variables such as 


a6 


Specific humidity and LCL, nite  @re"Luncti ons of air 
temperature. These variables are fundamental to the deter- 
Mination of the rate of change of all of the model parane- 


ters. 
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RScuEe i Mixed Layer Depth 
Be crved vs Model (Case [ 


Cbservations indicated a shallowing of the mixed layer 
twomtigcomenconigh 2300, Figs. 3.7. This was followed by a 
gradual deepening. It is due tc the surface heating and the 
distribution of heat by wind mixing followed by radiative 
cooling and convective overturning. The model predicted a 
continuous deepeniny of the mixed layer due to the substan- 
tially reduced values of shortwave radiation and the resul- 


tant sea surface cooling and convective instability. 


OFF 


D. CASE II: 1554 PDT O04 NOVEMBER-1600 05 NOVEMBER 1983 


The initializing variables for the second 24 hour period 
are listed in table II. Uncoupled and coupled model cutputs 
are presented in Figs. 3.8 and 3.9 respectively. 


}~ 


1. Uncoupled Mode 


SS = 


In this case, (Fig. 3.8), the cloud development is 
less extreme than that which arose in Case I. A cloud of 
minimum thickness formed when the model was initiated {20 
m) . As the run progressed, the model determined a fairly 
coms tant Lely whereas inversion height increased rapidly 
from 1600 through (900 00(492 te) 76o en The inversion 
attained a height of 1093 mat 0630 and remained at this 
Jevel until 1030, at which time it began to lower. At the 
conclusion of the model run the inversion had dropped 23 n 
from its maximum of 1093 ao. The LCL decreased slightly 
throughout the period. It began at a level of 472 fm and 
reached a final height of 422 np. The cloud thickness hada 


Maximum of 701 m at 0930, 05 November. 


The behavior of the inversion is reflected in the 
modeled entrainment rate. As the inversion height increased 
in the early portion of the _ period, entrainment rate 
attained an absolute ma Xi mum value of 3.68 cm/s. 


Thereafter, the predicted rate of entrainment decreased, 


ultimately attaining a value of .679 cm/s. The final value 
represented an absolute mininmun. Since the prescribed 
subsidence rate was held constant, the rate of entrainment 


again dominated the evolution of the inversion height. 
The LCL remained relatively constant throughout. 
This was a result of the compensating effects of dry air 


entrainment and an increasingly larger air-sea temperature 
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Figume 3.9 Coupled Model-Case II 
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TABLE II 


Initialization Yalues-Case fI 


Julian Date Sfc Pressure (mb) Sea Sfc Temp (°C} 
308 1G22:..3 16.70 

Inversion Ht (nm) Mixed Layer Potential Temp (°C}j 
G92 80 

Air Temp Jump (°C Temp Lapse Rate (°C/n 
4°EP p ( C) P wap se ( ) 


eee (g/kg) Spec ect ee Jump (u/kg n) 


Spec_ Humidity Lps Rate (g/kg n) ner 


0017 048) 12 
Subsidence (n/s) Mixed be ice Depth (m) 
0039 5 
Ocean eee Jump (°C) Ocean oP ee tet (°C/nm) 
Wind: Time Speed (n/s mnie © 
1600 P eu 340 >! 
18 30 6a5 340 
27100 5.6 Boul 
23.30 36 5 BUS 
0 200 Sue 325 
04 30 23 Q16 
07900 Dees SG 
0930 4.8 333 
1200 Doe 270 
1430 2.4 265 





dirierence. The former would force an increase in the LCL, 
the latter a decrease. This balance is also evident ir the 
nearly constant predicted specific humidity values. AS the 
moisture flux into the layer increased (due to the 
increasing air-sea temperature difference) the incorporation 
of dry, warm air (entrainment) resuited in a slight varia- 
tion iN speemenc humidity. 

Air temperature in the mixed layer, if viewed in 
EcullomOrperaclataive heating and cooling, woulda be expected to 


decrease after sunset and increase after sunrise. Ths. 1S 


4 4 


the case as shown in Fig. 3.8. The modelea shortwave radia- 
tion 1S again underestimated. As a result, a Significant 
decrease in air temperature occurred through 16 hours into 


the period due to cloud top long wave radiation. 


Z. Coupled Model 


As observed in Case I, the coupied and uncoupled 
model predictions are essentially identicai in this case. 
The coupled model prediction demonstrated a decrease in sea 
surface temperature and a deepening of the mixed layer. Sea 
surface temperature decreased from 18.7 to 18.63 °C. The 
moisture flux decrease, due to the reduced sea temperature, 
was too small to force any appreciable change in specific 
hums it y Ob cee Tae deepening of the nixed layer from an 
initial depth of 35 to a 37.4 m had no effect upon the 
atmospheric boundary layer. 

While Case I and Case II experienced similar 
decreases in sea surface temperature, the former deepened 
the mixed layer by a factor of two greater than the iatter. 
However, friction velocity was larger in Case I than Case 
Lie Therefore, the variation of the mixed layer depth was 


in response to wind forcing. 
3- Qbservations vs Model Predictions 


The predicted and observed specific humidity were 
significantly different, sFig aoe There is no difference 
between coupied and uncoupled predictions as they exhihit a 
Slight decrease throughout the period. Ovserved specific 
humidity varied in acyclical fashion due to the relative 
AUNLaaty varlarpon. 

Values of specific humidity were caiculated exter- 
hally using air temperature, pressure and relative humidity. 


Unlike Case I, a large decrease in relative humidity was 
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Bee 3610 Specific Humidity 
bserved vs Model (Case II) 


experienced. Specific and relative humidity changes are 
forced by surface moisture flux and entrainment. Since the 
predicted and observed values for air and sea sucface 
temperature are nearly the same, the difference lies in the 
entrainment. As noted earlier, the predicted rate of 
entrainment decreased through the period. AGfia. Cit barn— 
ment must have increased during the siddie of the run and 
decreased at the end. 

The radiosonde data (Fig. 3.11) support the observed 
decrease in mixed layer specific humidity. Mixeac layer 
specific humidity begins at a value of 9.76 g/kjy, decreased 


to 7.84 g/kg 16 hours later, and attained a value of 8.98 
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wegezs nNOours into the period. The predicted entrainment 


1s, therefore, smaller than that observed. 


LEGEND 
OBSERVED 
UNCOUPLED 
COUPLED 


= 
® 
O 
O) 
Q@) 
oo 
® 
a 
= 
ae, 
O 
cc 
@ 
2 
= 
® 
Es 
SS 
-) 
ie 


{ 
2200 0400 1000 


Time 





Figure 3. 12 Air Temperature 
Observed vs Model (Case IT) 


The air temperature was in close agreement for 
observed, uncoupled and coupled cases exists 16 hours into 
tne period, Fig. 3.12. From this time, the predicted values 
increase slightly while the observed aic temperature 
increases bv greater than 0.5°C. At the end of the period a 
temperature difference of 0.86°C exists. 

Predicted sea surface temperature agrees with 
observed values 20 hours into the run, Fig. enc Se The 


0 
okserved sea surface temperature increased 0.3 C from 16 to 
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Figure 3.13 Sea Surface Temperature 
Observed ws Model (Case ITI) 


24 hours after initialization (1200 to 1600 local time) due 
to shortwave radiative heating. Predictions fail to demon- 
strate this increase due to the extensive cloud cover and 
associated shortwave radiation deficit. This dichotomy is 
absent in the early stages of the run since the model was 
initialized with the observed sea surface temperature. 
Shortwave radiation 1S again underestimated, Fig. 
3.14. A clear sky was observed as depicted by the agreement 
of observed and clear sky values of shortwave radiation. 
The model predicted a thick cloud layer and assumed horizon 
to horizon coverage. This prediction resulted in a severe 


underestimation of shortwave radiation at the sea surface. 
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Figure 3.14 Shortwave Radiation 
Observed vs Model (Case IT) 


As in Case I, the effects of this radiation deficit are wide 


ranging throughout the model. 


Through the comparison cf coupled model and observed 
variations of the mixed layer depth, a periodic retreat of 
the mixed layer was observed, Fig. 3.15. A minimum depth 
was attained eight to ten hours after the maximum in short- 
wave radiation due to the delayed oceanic response to 
heating and the relatively light winds forcing the mixing. 
Additionally, the effects of tides and internal waves are 
most likely present as evidenced by the periodic variation 
in mixed layer depth. The shortwave radiation deficit 


forced a deepening of the mixed layer, which resulted from 
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Figure 3.15 Mixed Layer Depth 
Coupled vs Model 


cooling of the sea surface and convective mixing. Thies 
pattern corresponds to the remainder of the modeled 
variables. However, it is significantly different than the 


res;,onse of the ocean to the atmospheric forcing. 


E. CASE III: 1600 03 MOVENBER- 1600 05 NCY EMBER Mioes 


TO compare model performance and observations for a 
period greater than 24 hours, a 48 hour period was selected. 
Due to the limitations imposed upon the experiment by the 
eivironment and equipment malfunction, the time frames for 


Cases I and If were combined to provide the longer forecast 
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TABLE III 


Initialization Values-Case III 


Julian Date Ste Pressure {mb) Sea Ste Temp {~C) 
3085 1022.4 Lonwog 

Inversion Ht (nm) Mixed layer Potential Temp (°C) 
926 iy eat 

Air Temp Jump ('C Tonpebarseshateu, C/m 
10.3 peo 2 “Ob65 : 


eee eat Ot ty (g/kq) Spec Sane Jump (g/kg m) 


| 
Spec Humidity Lps Rate /kg a ine, (0 
Pp Se ng Y +p (g/kg @) noe ) 
Subsidence (m/s) Mixed Bier Deptn (m) 
= O52 U 
Ocean Zemp Jump (°C) Ocean Temp Gradient (°C/n) 
Wind: Time Speed (m/s) Dacwee 2) 
1600 5.4 340 
18 30 6.5 340 | 
2100 5210 a" | 
2330 5.8 345 
0200 Beye 225 | 
0430 eS 016 
0700 28 344 
09 30 eno oo | 
12100 Dye 2 270 
1430 ea 265 
[ 
period. The output variables from the last tine step of 


Case I were employed as initializing variables for the 
second 24 kour period of Case IIi. These variabies are 
lieted in-vable III. 


The most obvious difference in the Case III prediction 
as compared with that of Case II was that there were no 
rapid changes in the predicted variables, Figs. 3.16 and 
327s The most obvious difference He the SCase Tit 


PREe@lCtlOuyas Geompared With that of Case II was that there 
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Figures 3. 16 Case III Uncoupled Model 
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Figure 3.17 Case III Coupled Model 
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were no rapid changes in the output variabdies. The inver- 
sion height increased gradually from an initial value of $28 
m to a final value of 1146 on. The LCL ranged from anh 
initiai value of 464 m toa Miriaum value of 438 m toa 
Tinai height of 78eemen- The LCL and inversion height 
attained final values of 424 and 1073 m, respectively, inf 
Case II. Cloud thickness changed slightly during the 
period. It began at 464 m and increased to 663 mn. Cloud 
thickness for Case II grew to a thickness of 649 oa. The 
radiosonde soundings for the latter 24 hours period of Case 
Tie Pane 3.11) revealed the height of the inversion to be 
increasing. Though the profiles are not ideal, the final 
inversion base height is of the order of 1900 a. Hence, 
fairly yood agreement exists fer Cases II and III in this 
regard. 

Temperature also varied gradually during the period. 
Both Cases II and III demonstrated a diurnal variation. Air 
temperature profiles are depicted in Figs. 3.18 ani 3.19. 
Close agreement exists from the uncoupled and coupled 
models. There 1s essentially no difference between the 
output for these cases. The shortwave radiation deficit and 
associated lack of heating was obvious from 16 to 24 hours 
into the runs. An agreement of temperature values at the 
beginning of eacn run occurred as the result of initializa- 
tion and because the runs began at 1600. Darkness occurred 
Shortly thereafter. Thus there waS a Minimal amount of 
shortwave radiative heating. 

A very close agreement in specific humidity resulted for 
both uncoupled and coupled models, and for the 24 and 48 
hour periods, Figs. “3220, andeeoece AS highlighted in the 
discussion of Case II the difference between observed and 
modeled specific humidities was due to an apparent increase 
in the actual rate of entrainment of warm dry air into the 
mixed layer, while modeled entrainment remained fairly 


constant and smaller in magnitude. 
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Figure 32°76 Air Temperature Comparison 
Uncoupled Models 
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Figure 3.20 Specific Humidity 
Uncoupled Model 
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Shortwave radiation varied as expected for the 48 hour 
Boe owweancluding the@latter 24 hours, Fiys. 3.22 and 3.23. 
The order o£ magnitude of the difference between observed 


and modelled values has been the same for the three cases 
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analyzed. As in Cases I and II, the extreme cloud cover and 
Minimal shortwave radiation have forced significant differ- 
ences in modeled variables. The impact of the predicted 
radiation deficit was obvious in the behavior of air temper- 
ature in Case III and was also evident in the difference 
between observed and modeled mixed layer depth and _ sea 
surface temperature. 

Rhile the mixed layer deepens in Cases II and III, the 
observed mixed layer depth resfonds as expected to surface 
heating and wind mixing, Fig. 3.24. As described in Cases I 
and II, an eight hour lag exists between maximum shortwave 


radiation and maximum retreat of the mixed layer, which 
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Figure 3.23 Shortwave Radiation 
Coupled Model 


results from the relatively slcw response of the ocean to 
atmospheric thermal forcing. 

The modeled cooling of the sea surface temperature in 
all cases is yet another effect of the shortwave radiation 
deficit, \Fageecem>. 


Fe. EXPERIMENTAL SIMULATION 


Based upon the results of this study, research into the 
problem of cloud thickness prediction and shortwave radia- 
tion deficit has begun. Case I variables (Table I) were 
used to initialize a revised version of the model, Figs 3.26 
and. 27 Of major concern is the difference in the 
predicted evolution of cloud thickness. The uncoupled and 


coupled models predicted clouds of thicknesses ranging from 
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Pigure 3.24 Mixed Layer Depth 
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Figure 3.25 Sea Surface Temperature 
Coupled Model 
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FaGcune 3.27 Coupled Model Revised 
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50 to 100 a. Case™i clouds were “orlyinally Vprediee ie 
attain a thickness Oiieeein (tec eo) Predicted values 
f shortwave radiation are in better agreement with observa- 
tions (aS compared with prior cases) although differences do 
exist, Fig. 3.28. Mixed layer depth was predicted to shal- 
low whereas the earlier cases demonstrated only deepening. 
This too was a result of the difference in the predicted 


evolution Ot ShHOEUWave LTadzat2on- 





pa 600.0 
£ 550.0 LEGEND 
9 Bea B @lsisie a0 
a + Cale tee) 
2 450.0 x COUPLED 
~ 400.0 
5 
me S50e8 
Gu 30080 
) N 
OU 25046 
O 
> one | if 
qa) 
> 15028 
5 
as 100.0 - tp 
a 
O 5056 av 7 ad 
Re 
i) G18 barre Orin FITIVIYITITOTUT ITI YITIYIYIT.1 “4 
1600 2200 O400 1000 1600 
Time 
FL0UreeS3e7o Shortwave Radiation Comparison 


The evolution of predicted inversion hergnc. LGl eae 
temperature and specific humidity by the uncoupied snodel 
demonstrated an extreme jump 20 hours after initialization 


{1200 04 Nov). AS demonstrated previously, dynamic feedback 


6 0 


continued to occur anong the predicted variables. 
Variations of such magnitudes however, over a period of less 


than 30 mirutes are unrealistic. 


6 1 


Uncoupled and coupled model predictions overestimate 
cloud thickness and assume an unrealistic horizon to horizon 
cover. AS a result, short and long wave radiative effects 
are minimized. The impact of this discrepancy affects all 
of the model's physics as the fluxes of heat and moisture, 
variation of sea surface temperature and mixed layer depth 
are improperly forced. Hence, the values for boundary layer 
variables such as specific humidity, air temperature, LCL 
and inversion height differ fron observations. 

Given this shortcoming, the model demonstrates little 
variation between the coupled and uncoupled tnodels. The 
value of oceanic-atmospheric coupling has been documented 
(O'Laugian 332). The model behaves well for an extended 
period (48 hours) when reinitialized after 24 hours. 

The preliminary results obtained through the use of a 
modified version of the model indicated that the model will 
predict commonly observed variations in parameters such as 
mixed layer depth. Similar predictions can be obtained 
through the use of different expferinental initializing vari- 
ables. The key question which must now be investigated is 
why did the original coupied and uncoupled models behave as 
they did with the MILDEX observations. Due to the tentative 
nature of these results, definitive conclusions concerning 
the radiation and cloud modeling as it now exists cannot be 
drawn. The MILDEX data set 1S somewhat limited in its 
applicability cue to gaps in the many variables reguired to 
initialize the model. The main advantage of the data set 
was the lack of platform motion during the exercise, thus 
limiting advective effects which such motion may introduce. 
The data set should continue to be valuable in future 


studies with this model. 
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Cloud physics medeling must continue to Le reviewed 
along with verification and sensitivity analyses. eee 
aayer cloud model should be incorporated to examine radia- 
tion prediction. Though adeéguate, the slab approach 
currently employed is no longer essential as a means of 
reducing computation time. Adit rona ll ys a Gavpabilitywete 
account for partial cloudiness is absolutely essentiai. The 
total cloud cover assumption currently employed nas been 
demonstrated by this study to be the most serious short- 
coming of an otherwise rundanentally sound model. A statis- 
tical approach in the description of cioud cover percentage 
is suggested as a pessible solution. Further study is 
essential in this area. 

The model should remain in a research phase of deévelop- 
ment and not be distributed for operational evaluation until 
the probiems outlined herein are rectified. Once the model 
has peen further refined to alleviate these probiems, it 
Shouuad tnen be incorporated into a tactical environmental 
prediction model. Then the model will be a useful tool for 
the Gceanography Officer in supfort of fleet operations. 

As for data acquisition, the performance of a navaid 
radiosonde system should be analyzed. The inability to 
record pressure and altitude for several launches severly 
restricted the amount of useful radiosonde data in MILDEX. 
Additionally, graduate student involvement in data acquisi- 
tion e€xperiments is essential. Such an individual would be 
able to acquire a background on the data base and the atmos- 
pheric and oceanic Situation. An exercise conducted in an 
area such as Monterey Bay would be useful. Such a limited 
cody of water would further limit the advective effects of 


the ofen ocean. 
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